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We theoretically study the decoherence and the dynamical decoupling control of nitrogen-vacancy center 
electron spins in high-purity diamond, where the hyperfine interaction with ^^C nuclear spins is the dominating 
decoherence mechanism. The decoherence is formulated as the entanglement between the electron spin and 
the nuclear spins, which is induced by nuclear spin bath evolution conditioned on the electron spin state. The 
nuclear spin bath evolution is driven by elementary processes such as single spin precession and pairwise flip- 
flops. The importance of difl'erent elementary processes in the decoherence depends on the strength of the 
external magnetic field. 
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I. INTRODUCTION 

The optically detected magnetic resonance of single 
nitrogen- vacancy (NV) centers in diamond^ triggered the re- 
search on spin coherence of single NV centers for applica- 
tions in quantum information processinj^and ultra- sensitive 
metrology.^ As one of the most important advantages, NV 
center electron spins have long coherence time up to the mil- 
lisecond 

timescal^^Stll 

even at room temperature. The long 
coherence time in such systems makes them an ideal platform 
for studying quantum sciences and technologies. Understand- 
ing the decoherence mechanisms and prolonging the coher- 
ence time of NV center spins are of fundamental i mporta nce 
in exploring a number of surprising physical effectPm and 
novel applications .^^^[^ 

In general, electron spin decoherence in solids can result 
from phonon scattering through spin-orbit coupling, interac- 
tion with paramagnetic impurities, and hyperfine interaction 
with nuclear spins. Being a light-element material, diamond 
has weak spin-orbit interaction, and therefore the spin-lattice 
relaxation of NV center spins caused by phonon scattering 
has rather long timescales, Ti =7.7 ms at room temperature 
and longer than seconds at low temperature.^^ In most cases, 
the NV center electron spin decoherence is not limited by the 
spin-lattice relaxation. The coupling to electron or nuclear 
spins in diamond can induce faster decoherence. In type-lb 
diamond samples, the main paramagnetic centers coexisting 
with the NV center are nitrogen donors with one unpaired 
electron spin (the PI centers). These paramagnetic centers 
couple to the NV center electron spins through the dipolar 
interaction, and form an electron spin bath. The NV center 
electron spin coherence time is inversely proportional to the 
concentration of the PI centers. For a PI concentr ation 
-10^ ppm, the NV center spin coherence time is ~ jds^^ 
In Fig.[T] we plot the estimated NV center free-induction de- 
cay (FID) time r*^^ as a function of the PI concentration, 
which is the lower iDOund decoherence time contributed by the 
electron spin bath. When the PI concentration is decreased 
to < 10 ppb, as in the high-purity type-IIa samples, the decay 
time T* due to the electron spin bath will exceed millisec- 
ond. In this case, for NV centers in diamond with natural 



abundance of ^^C, the decoherence will be dominated by the 
hyperfine interaction with the ^^C nuclear spins, which form a 
nuclear spin bath [see Fig.[2jb)]. In this paper, we will focus 
on the NV center electron spin decoherence in the ^^C nuclear 
spin bath. 

In order to protect the electron spin coherence from the en- 
vironmental noises, dynamical decoupli ng (pD)^" ^ control 
has been demonstrated as an eflftcient way! 1 1 ^ ^^ In general, 
under DD control, the electron spins are repeatedly flipped, so 
that the eflTect of the noises is averaged out, and the coherence 
time is prolonged. In this paper, the NV center spin coher- 
ence under DD control is investigated, and the basic physical 
processes in the nuclear spin bath are analyzed to reveal the 
decoherence mechanisms in diflTerent magnetic field regimes. 

Similar problems of electron spin decoherence in nuclear 
spin baths have been investigated in other systems such 
as electron spins in quantum dots (QDsj^^^ and e lectron 
spins of shallow donors in silicon, such as Si:F *^^l^^^ and 
Si:Bi!ailThe Spin decoherence in diamond has several dis- 
tinct features as compared with that in QDs and shallow 
donors. 

First, NV centers are deep-level defects in diamond. In QDs 
and shallow donors, the electron wave functions extend to a 
few or even tens of nanometers, and the couplings between the 
electron spins and nuclear spins are dominated by the isotropic 
Fermi contact interaction. In NV centers, the electron wave 
function is localized around the vacancy site, extending only 
several angstroms. The Fermi contact part vanishes quickly 
as the distance increases. Also, the natural abundance of the 
^^C, the only C isotope with non-zero spin, is only about 1.1%. 
The electron spin couples to the bath ^^C spins mainly through 
the dipole-dipole interaction. The Rabi oscillation of central 
spins in dipolar-coupled spin baths was studied, and complex 
decay behavior of central spins under continuous driving was 
revealed in Ref. ^ Here we will focus on the spin deco- 
herence behavior under pulsed-DD control. The anisotropic 
nature of the hyperfine interaction plays an important role in 
the NV center decoherence. 

Second, the bath spins involved in the decoherence of NV 
centers are much fewer than those in QDs and shallow donors. 
In QDs and shallow donors, a large number of nuclear spins 
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FIG. 1. NV center electron spin decoherence in an electron spin 
bath. With only the electron spin bath considered, the FID decay 
time T* ^ ^ is inversely proportional to the concentration of the bath 
electron spins, and exceeds ~ 1 ms for the concentration lower than 
0.01 ppm (the shadowed region). 



(e.g., ~ 10^ for Si:P with natural abundance ^^Si, and 10^ ~ 
10^ for QDs depending on the QD size) contribute to the deco- 
herence. The contribution of each nuclear spin (or spin pair) 
to the total decoherence is small and can be well described by 
the nuclear spin dynamics in the short time limit. In contrast, 
the nuclear spin baths of NV cente rs in diamond consist only 
hundreds of ^^C nuclear spin^^, which gives rise to a long 
coherence time violating the short time condition for the nu- 
clear spin dynamics. This brings th e uni que features in the 
NV center electron spin decoherence ] 

Particularly, due to the random location of the ^^C nuclear 
spins on the diamond lattice sites, the nuclear spin baths of 
NV centers could be inhomogeneous, i.e., the distance from 
a certain ^^C spin to its neighboring spin can have a large 
fluctuation around the average distance between neighboring 
spins. This distance fluctuation enables closely bonded nu- 
clear spin clusters to appear around the NV centers. Due 
to the inverse-cubic dependence of the dipolar interaction on 
the inter-nucleus distance, such closely bonded clusters could 
have much faster dynamics than those weakly bonded nuclear 
spins, and therefore be distinguished in the decoherence pro- 
cess. As an example, in Ref. [TT] we have demonstrated that 
a nuclear spin dimer, i.e. two nuclear spins occupying a C-C 
bond, can induce characteristic oscillations on the decoher- 
ence profile. Here, we will further show that the quantum dy- 
namics of closely bonded nuclear spin clusters [see Fig.[2jb)], 
as localized elementary excitations in the interacting nuclear 
spin bath, will become the main deocoherence mechanism un- 
der certain magnetic fields. 

Theoretically, to investigate the electron spin decoherence 
in an interacting nuclear spin bath, one needs to solve the 
quantum many-body dynamics of the spin bath. Various 
methods have been developed to treat this kind of prob- 
lem, including the cluster expansion,^ the pair-correlation 
approximation,'^^'^^' and the disjoint clusters approach. 
When investigating electron spin coherence of NV centers 
under DD controls, one usually has to take into account the 
higher order correlations and the dynamics beyond the short 
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FIG. 2. (a) Schematic of the NV center structure. The magnetic 
field is assumed to be applied along the direction pointing from the 
nitrogen to the vacancy (the [111] direction), (b) Schematic of the 
NV center in a ^^C nuclear spin bath. The ^^C atoms with nuclear 
spins (the red arrows) are randomly located on the diamond lattice 
(the gray circles). A nuclear spin dimer is encircled, (c) Schematic 
of the DD control sequence. The coherence is generated by the initial 
nil pulse (red squares). The electron spin is repeatedly flipped at 
times tk by the 7i pulses (blue squares). Finally, the coherence is 
converted to the level population by the last Jijl pulse. 



time limit. I n thi s case, the cluster-correlation expansion 
(CCE) methoc?^^ previously developed in our group is es- 
pecially suitable, which has been proved to be an efficient 
method in treating the decoherence in correlated nuclear spin 
baths of finite size.^^ This paper focuses on the decoherence 
effects of different elementary processes in the nuclear spin 
baths under various conditions. Comparison of different theo- 
retical methods is beyond the scope of this paper. 

The paper is organized in the following way. Section |II| 
gives a description of the NV center electron spin decoher- 
ence problem and the CCE method. 



we an- 



In Section III 

alyze the basic spin processes involved in the decoherence. 
The results of NV center electron spin coherence in different 
magnetic regimes under DD controls are discussed in Section 
IV Section [V] gives the conclusion. 



II. MODEL AND METHOD 



A. System Hamiltonian 



In this paper, we consider negatively charged NV centers 
in diamond. The electronic ground state of the NV center is 
a triplet state with the spin quantum number 5" = 1 . The NV 
center spin is coupled to the ^^C nuclear spins {IJ, which are 
spin- 1/2 's with the natural abundance /?nat = 1.1%, randomly 
distributed on the diamond lattice. The full Hamiltonian of 
the whole system is written as 



H - //mv + ^bath + fi\ 



(1) 
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The Hamiltonians of the NV center and the bath spins in mag- 
netic field B are 



ifNV = -reB•S + A^^^ 

^bath = -TnB • ^ 1/ + if dip, 



(2) 

(3) 



where = -1.76 x 10^^ rad s"^ T"^ and yn = 6.73 x 
10^ rad s~^ T"^ are the gyromagnetic ratios of the electron and 
^^C nuclear spins, respectively, A = 2.87 GHz is the zero-field 
splitting of the electron spin, and the z-axis is taken along 
the direction pointing from the nitrogen to the vacancy [see 
Fig. [2ja)] unless specified otherwise. In this paper, the mag- 
netic field direction is assumed along the z direction. We fo- 
cus on the dependence of the decoherence behavior on the 
magnetic field strength. The eff'ect of the magnetic field direc- 
tion on the electron spin coherence time has been studied in 
Refs.[ni[38land|4ll 

In Eq. ([3]), the Hamiltonian //dip includes the dipolar inter- 
actions between nuclear spins 



//c 



dip 



3(lrr,;)(r,,.-I,.) 



(4) 



where r^y is the displacement from the ith nuclear spin to the 
jth nuclear spin, and the interaction strength is characterized 
by Dij = [loy^l (Anr^^j) with yUo being the vacuum permeability. 

The electron and bath spins are coupled through the hyper- 
fine interaction, which is described by 



/i,nt = S-^ArI/, 



(5) 



where A/ is the hyperfine interaction tensor for the /th nu- 
clear spin, including the isotropic Fermi contact part and the 
anisotropic dipolar interaction part. The Fermi contact part is 
important for nuclear spins near the electron spin. The ^^C 
nuclear spins appearing in the first few coordinate shells will 
induce fast (> MHz) electron spin echo envelope modula- 
tion (ESEEM)."^^ Instead of the ESEEM, this paper focuses 
on the overall decoherence, which arises from the coupling to 
the large number of nuclear spins relatively far away from the 
NV center. For these nuclear spins, the Fermi contact is not 
important, and the hyperfine interaction takes the dipolar form 



An 



(6) 



where r/v is the displacement of the /th ^^C spin from the va- 
cancy site. 

B. Electron spin coherence 

The total Hamiltonian in Eq. ([T]) is rewritten as 

// = //Nv + b-S + //bath, (7a) 
b^^A,-.I,-, (7b) 



where the influence of the hyperfine interaction in Eq. ([5]) is 
expressed in terms of an efl'ective field b, provided by the nu- 
clear spins, coupling to the NV center spin. Either the thermal 
distribution of the nuclear spin states or the dynamics of nu- 
clear spins causes the fluctuation of the efl'ective field. Thus, 
the eff'ective field b is regarded as the noise field, which in- 
duces the decoherence of the NV center spin. 

For the magnetic field applied along the z direction, the 
magnetic quantum number m is a good quantum number, and 
the eigenstates of //nv are denoted as \m) with m = 0, and 
±1, and the corresponding eigenenergies ojm = m^A - my^B. 
With these eigenstates, the hyperfine interaction in Eq. ([5]) is 
expanded as 



+1 



^int= Yj I^><^I^^'^'«' (^^) 
m,n=-\ 

bm,n — ^m,n ' ^ = ^m,n ' A/ • I/, (8b) 



where Sm,n = <m|S|^). The noise operator bm,n contains only 
the nuclear spin operators. 

For the NV center spin, the diagonal part of the noise oper- 
ator bm,m induces the energy shift of the level m, and the off'- 
diagonal part bm,n formi^n causes the transition between |m) 
and 1^). Since the NV center has a zero-field splitting in the 
order of GHz, which is much larger than the typical hyperfine 
interaction strength (< MHz), the electron spin can hardly be 
flipped by the hyperfine interaction. As a higher-order eff'ect, 
the hyperfine interaction can still induce indirect coupling be- 
tween nuclear spins by virtual flips of the electron spin. This 
indirect coupling can afl'ect the electron spin coherence only 
when the nuclear spins involved in the virtualprocesses are 
strongly coupled (> MHz) to the electron spin.^Our previous 
stud>^^ has shown that, in most cases that the NV center elec- 
tron spin levels are far from the degenerate point, the indirect 
coupling is much smaller than the intrinsic dipolar interaction 
in Eq. Q, and has negligible efl'ect on the NV center spin 
decoherence. As a consequence, the ofl'-diagonal terms with 
m n can be safely neglected, and a pure dephasing Hamil- 
tonian is obtained as 

+1 

H^H^y^ b,S, + //bath = Yj I^X^I ^ 



(9b) 



where H^'^^ is the bath Hamiltonian conditioned on the elec- 
tron spin state |m). The operator bm,m is written in the follow- 
ing form 

bm,m = ^ m (z . A,-) • 1/ = ^ Af'^ . I,-, (10) 

where A^"^^ = m (z • A/) is the efl'ective hyperfine field for the 
ith nuclear spin when the electron spin is in the state |m). 



The hyperfine field in Eq. (10) contains two parts, namely, 
the isotropic and anisotropic interactions, according to 
whether or not the nuclear spin quantum number is conserved. 
The isotropic part contains the coupling along the z direction 



4 



(i.e. the A^^hf term), and provides frequency shifts to nu- 
clear spins. The anisotropic part contains the coupHng along 



the X and y directions (i.e. the Af^Jlf and A^^^f^ terms), and 
causes the nuclear spin flipping. Since the nuclear spin flip- 
ping process does not conserve the Zeeman energy, the efl'ect 
of the anisotropic coupling on the decoherence depends on 
the magnetic field strength. The role of anisotropic coupling 
under diff'erent magnetic fields will be discussed in detail in 
SectionHniand Section Hvl 

In the absence of electron spin flipping process in the 
Hamiltonian Eq. ([9]), the population of each electron spin 
states will not change. The coherence of the electron spin 
at time t is defined by the average value of the transverse spin 
component as 



L(t) = 



(11) 



where 5^ = S x iSy and p(t) is the density matrix of the 
total system of electron and bath spins at time t. For the initial 
time ^ = 0, the system is prepared in a product state with 
p(0) = Pb ^ \i//q(0)}{i//q(0)\, where pe is the density matrix of 
bath spins, and the electron spin is in a superposition state of 
eigenstates \m) and \n), i.e. \i/^q(0)) = a\m) + p\n). 

For the temperature much higher than the nuclear spin Zee- 
man energy (<MHz ~ yuK), the density matrix pb of the bath 
spins is well approximated as pb =1/2^, where M is the num- 
ber of ^^C spins, and I is a 2^ x 2^ identity matrix. With N 
DD pulses applied at ^i, ^2, • • • , ^a^, which flip the electron spin, 
the coherence between \m) and \n) states is expressed by^ 



L(0 = TrB[---^-^^'''"^^ - -PB^" -e 



(12) 



where Tk = tk - h-i for = and ^a^+i = ^ is the ^th free 
evolution interval between the {k - l)th and ^th pulses [see 

In this paper, the pulses are assumed to be ideal, which flip 
the electron spin instantaneously with no error. We will focus 
on the periodic DD (PDD) pulse sequences, in which the elec- 
tron is periodically flipped at times tk = (2k - l)t/2N. This 
control sequence was shown to be an efficient way for pro- 
tecting the electron spin coherence of NV centers. The 
comparison and analysis of the eflSciency of various kinds of 
DD control schemes are beyond the scope of this paper and 
will be discussed elsewhere. 



C. Cluster-correlation expansion method 

The decoherence of an electron spin in an interacting bath 
is a many-body problem. The CCE method is employed to 
solve the problem. 

The key idea is that the decoherence function L(t) can be 
expressed as the product of cluster correlations 



Lit) 



Y]Lc(t), 
c 



(13) 



with the unfactorizable correlation Lc (t) of cluster C recur- 
sively defined by 



Lc(t) = 



Lc (0 



Uc'ccLc (t) 



where Lc (t) is calculated in the similar way to Eq. ( 12 ) 



Lc(t) = TrB 



(14) 



(15) 



where only the interaction within the cluster C is included in 
the bath conditional Hamiltonian H^^\ In realistic calcula- 
tion, the expansion is truncated at a certain size K of clusters 
(defined as the number of spins in the cluster), i.e. |C| < K, 
and the result is denoted as CCE-^. The final electron spin 
coherence is obtained by increasing the value of K until the 
results get converged. The CCE convergence order depends 
on the underlying microscopic decoherence mechanisms in 
the timescale of interest. Thus, CCE is not only an efficient 
calculation method, but also provides a tool for analyzing the 
dominating decoherence mechanisms. Typically, the trunca- 
tion up to four- spin clusters, i.e. CCE-4, can give a convergent 
result. For the electron spin coherence under many -pulse DD 
controls, as the coherence time is prolonged to ~ 10 ms, the 
calculations are convergent at CCE-6. 



III. PHYSICAL PROCESSES 

For NV centers in nuclear spi n bath s, the noises can be of 
either thermal or quantum nature .1^1^ 

At room temperature, the random orientations of the nu- 
clear spins result in thermal noises. The density matrix of 
nuclear spin is written as 



Pb 



1 

2^ 



(16) 



where |/) is chosen as the eigenstate of the noise operator 
[see Eqs. ^ and ([9])] with 



b,\J} = bj\J). 



(17) 



The thermal distribution of the nuclear spin states |/) leads to 
the random Overhauser field b j for the NV center spin, which 
is the thermal noise. The thermal noises are so strong that they 
induce the NV center decoherence in several yus. However, the 
thermal noises are quasi-static and can be completely elimi- 
nated by spin echo. Then, the dynamical noises caused by the 
quantum evolution of the nuclear spins become important for 
NV centers under DD controls. 

The dynamical noises arise from the quantum evolution of 
the nuclear spins. Since the noise operator b^ does not com- 
mute with the bath Hamiltonian //bath, i-e. [b^, //bath] ^ 0, the 
eigenstates |/) of the noise field are in general not the eigen- 
states of the total Hamiltonian. During the evolution of the 
system, the nuclear spins initially in state |/) will be brought 
to the superposition of eigenstates of b^ corresponding to dif- 
ferent Overhauser field bj. Thus, the evolution of the nuclear 
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(a) 




FIG. 3. (a) Effective fields hf^ for the jth single nuclear spin in a 

magnetic field B and a hyperfine field A^f"^ conditioned on the elec- 
tron spin state |m>. (b) Bloch sphere representation of the bifurcated 
evolution of the jth nuclear spin under a single pulse control (Hahn 
echo). The single nuclear spin is driven by the eff'ective fields h^^^ 

and h^^^\ The distance Sj between the trajectories measures the de- 
coherence due to the jth nuclear spin. 



spins causes the dynamical quantum noises of the NV center 
spin levels. In order to investigate the decoherence mecha- 
nisms under DD controls, it is critical to understand the el- 
ementary processes in the nuclear spin evolution. Since the 
hyperfine interaction and the nuclear spin dipolar interaction 
involve one or two nuclear spins, respectively, the basic phys- 
ical processes occurring in the nuclear spin bath are the single 
nuclear spin precession and the nuclear spin pair flip-flops. 



A. Single nuclear spin dynamics 



m, = or ±1, and m ^ m') gives the bifurcated paths of the 
nuclear spins, and the electron spin coherence Lj(t) is deter- 
mined by the distance 



(19) 



between the trajectories on the Bloch sphere,^^ where 
\i//'^^\t)} and \i//'^^'\t)} are the nuclear spin state driven by the 



eflfective field h^^^ and h^^ \ respectively. 



According to Eq. ( p3] ), the 7th single nuclear spin contribu- 
tion to the FID and Hahn echo of the coherence between |0) 
and I + 1 ) states is calculated aJ^ 



q(0) 



Lj,FiD (0 = COS COS + sin sin cos (fj. 



q(0) 



7,Hahn 



(t = 2t) 



,(+1) 
2 "J 



1-2 sin^ (fi j sin^ sm 

^^2 2 



(20a) 
(20b) 



where 6^j^^ = y^h^J^^r and (pj is the angle between h|^^^ and 
(see Fig.B. Under A^-pulse FDD (PDD-A^) controls, the 



(0) 



decoherence at ^ = 2Nt due to 7th nuclear spin is expressed 
as 



|1 -2sin^Qfsin^(^|, 



N = 2k 



U/,Hahn(2T)L,;2;^(4^T) + Leorr(T), N = 2k + 1 



(21) 



fork= 1,2,.... The angles a and 6 are defined as 



In this subsection, we investigate the quantum evolution of 
the single- spin clusters. The conditional Hamiltonian of the 
7th nuclear spin is 



J J 



I, 



(m) 



B 



■Af/rn. 



(18a) 
(18b) 



The geometric picture of the eflTective field h^^^ is shown in 
Fig. [Sj^a), and the Bloch sphere representation of the single 
spin precession is illustrated in Fig.jSj^b). The north and south 
poles represent the basis states | t) and \ i), respectively. For 
the electron spin in state |0), the hyperfine field vanishes, and 
the nuclear spin rotates about the applied magnetic field, i.e. 
h^.^^ = B. For the electron spin state | ± 1), the eflfective field is 



J 

,(±1) 



B 



■Ar>/r„. 



Because of the anisotropic nature of the hyperfine interac- 
tion, the eflTective fields h^^^ and h^~^^ are in general not paral- 
lel to each other. This is dififerent from the cases of the QDs 
and shallow donors. In those systems, the hyperfine inter- 
action is dominated by the isotropic Fermi contact part, so 
the diagonal hyperfine field is always in the same direction 
as the applied magnetic field. As a result, in QDs and shal- 
low donors, the single nuclear spin dynamics does not signif- 
icantly contribute to the electron spin decoherence. 

In the NV center system, the precession of the 7th nuclear 
spin about the unparalleled eflTective fields h^^^ and h^^ ^ (with 



tana 



cos 



(22a) 



cos^^^^ cos^^.^^^ ■ 



■sin6>f sin6/^^^^^cos^y. (22b) 



In Eq. ( [2T] ), the coherence under FDD-(2k +1) control is writ- 
ten as the product of Hahn echo and PDD-2^ control with an 
additional correction term Lcorr(T). The expression of Lcorr(T), 



the derivation of Eqs. ([2Q|)-([22|) and the geometric pictures can 
be found in Appendix. 

The single nuclear spin contribution to the electron spin de- 
coherence depends on the magnitude of the applied magnetic 
field. In the absence of magnetic field, the eflTective fields are 



= and h^.^^^ = -AS^'Vrn, and all the evolution opera- 
tors in Eq. ([15]) mutually commute. Therefore the single spin 
dynamics does not contribute to the electron spin decoherence 
at the echo time. For weak magnetic field (B < I Gauss), the 
precession angle 6^^ about the magnetic field is small within 
the relevant decoherence timescale (r <^ ms). In this case, the 
single spin dynamics of each nuclear spin gives a minor con- 
tribution to the electron spin decoherence, which is character- 
ized by the distinguishability between the two trajectories, 
i.e. S^jOcylB^T^ <^ 1. 

In the opposite limit that the nuclear spin Zeeman energy 
much is larger than the hyperfine interaction, the single nu- 
clear spin contribution to the decoherence is suppressed. Ac- 
cording to the Bloch sphere picture shown in Fig. [3jb) and 
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FIG. 4. (a) Energy levels and eigenstates of two-nuclear spin Hamil- 
tonian H^^^.^ under zero magnetic field for the NV center spin in state 
|0>. (b) The energy shift and the transitions induced by the hyperfine 
fields. The transitions induced by the average hyperfine field A*^-^^ 
(the red thick arrows) are much stronger than those induced by the 
hyperfine field diff'erence a*^-^^ (the blue dashed arrows). 



1.0 




O 

c 

CD 
&_ 



_c 
o 
O 



0.5 



0.0 



(a) 



1.0 



CD 
O 

c 

CD 
i_ 

CD 
O 

O 



0.5 



0.0 



(b) 



PDD-50 
dinner 

exact 

Q approx. 



PDD-50 
separated pair 

exact 

Q approx. 




10 20 30 

Time (ms) 



Eqs. (20)-(22), the single spin contribution to the decoher- 
ence is controlled by the angle (pj between the two eflTective 
fields. Increasing the magnetic field so that y^B ^ "^j^^^' 
pseudo-fields are approximately parallel with each other, i.e. 
sin ifj ^ 1 . In this case, the two bifurcated evolution paths 
do not have significant distinguishability. In other words, the 
decoherence due to single nuclear spin dynamics is small in 
the strong field limit. 

For a medium magnetic field, the single nuclear spin pre- 
cession induces the oscillations of the electron spin coherence. 



As shown in Eqs. (20)-(|22]), there are two recovery periods of 



the total evolution time Ty = ANnliyJi^J^"^) for m = and 
+ 1. These two recovery periods correspond to the nuclear 
spin precessing about the eflTective fields h^^^ for full circles. 

Since the effective field h^^^ = B is actually the applied mag- 
netic field, the recover period tJ^ = 4Nn/(y^B) is the same 
for all the nuclear spins. The simultaneous recovery of coher- 
ence for all nuclear spins results in the periodical coherence 
revival eflTect, which has been experimentally observed in the 
Hahn echo case.l^^^ 



B. Nuclear spin pair dynamics 



In a two-spin cluster, the two nuclear spins are correlated 
due to the dipolar interaction. In general, there are two types 
of processes driven by the dipolar interaction, namely, (i) the 
non-secular spin flipping [e.g., the processes described by 
terms like /^/| and I^I^], and (ii) the secular spin flip-flop 
[e.g., the process described by ^^^J]. The non-secular flip- 
ping does not conserve the nuclear Zeeman energy while the 
secular flip-flop does. Consequently, they have difl'erent dy- 
namic behaviors under difl'erent magnetic fields. 



FIG. 5. (a) The electron spin decoherence between | + 1> and | - 1) 
states induced by a nuclear spin dimer under PDD-50 control and 
zero magnetic field. The black line is the exact calculation result us- 
ing Eq. (T5| l, and the red line with symbols is obtained by the pseudo- 
spin model [Eq.([26|]. (b) The same as (a), but for the two nuclear 
spins separated by 2.52 A (the next nearest neighboring pair). Notice 
the different timescales of the horizontal axes. The separated pair has 
much slower dynamics than the dimer. 



1. Non-secular spin flipping 

The non-secular process occurs when the nuclear Zeeman 
energy cost is less than or comparable to the dipolar interac- 
tion strength, i.e. y^^B < Dfj. Since the dipolar interactions 
between ^^C nuclear spins are at most ~ kHz (for the closely 
bonded pairs), the non- secular process will be activated in the 
weak magnetic field < 1 Gauss. 

For the sake of clarity, we consider first the zero magnetic 
field case. For the electron spin in the |0) state, the Hamilto- 
nian of the nuclear spin pair only contains the dipolar interac- 
tion. By choosing the quantization axis along the direction of 
relative displacement of the two spins, i.e. z || r^y, the Hamil- 
tonian H^^\, is written as 



(0) 
ij} 



Di 



[(7-7J+7r/;)-4/^/]], (23) 



where I- = R ± if. . Since the Hamiltonian H'f^\, commutes 
I I I {i,j] 

with the total spin angular momentum (1/ -h Ijf- of the two 
^^C nuclear spin-l/2's, the eigenstates of are classified 

as the triplet states, = \ TT>, Fo) = (I M) + I Ti)) / V2 
and = I U), and singlet state |^) = (| Ti) - I Ti» / V2, 
according to their total angular momentum. The eigenstates 
and the corresponding energy levels are shown in Fig.|4ja). 

For the electron spin in the | ± 1) states, besides the dipolar 
interaction between nuclear spins, the hyperfine interactions 
provide eff'ective fields for the nuclear spins. The conditional 
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(a) 



^ ^ decoupled ^'^^ 



pseudo-spin 




J_L i 



secular flip-flop 



decoupled 




FIG. 6. (a) Energy levels of two-spin cluster under a strong magnetic 
field. The two polarized states are decoupled and the dynamics of 
the two unpolarized states is described by a pseudo-spin, (b) The 
bifurcated evolution paths of the pseudo-spin on the Bloch sphere. 



The pseudo-spin is driven by the pseudo-fields and h^.^ 



(+1) 



The 



distance 6jk between the trajectories measures the decoherence due 
to the nuclear spin pair. 



Hamiltonian reads 



(l,-I,).(24b) 



In Eq. p4b| ), we have expressed the hyperfine interaction in 



terms of the average hyperfine field AjJ^jj = (A^^^^ -h A|^ 



and the hyperfine field diff'erence a^n = " ^ - 

In general, due to the presence of the hyperfine fields, the 
states l^") and |ro,±i) are not the eigenstates of H^fjy Because 
the coupling to the average hyperfine field A^f V conserves the 



total angular momentum, A|*jj^ will induce the transitions and 
the level shift within the triplet subspace. The coupling to 
the hyperfine difference a|*jj^ breaks the conservation of to- 
tal angular momentum, and causes the transitions between the 
singlet and triplet states [see Fig.|4jb)]. In particular, the tran- 
sition between l^") and l^o) states, which conserves the z com- 
ponent of the total angular momentum, is due to the secular 
spin flip-flop. 

For relatively near-neighboring nuclear spin pairs (with the 
inter-nucleus distance of several A), the average hyperfine 
field is usually much larger than the hyperfine field diflTer- 



ence, i.e. \^\'^^} \ » 



With this observation, the singlet 
state l^") is eflfectively decoupled from the triplet states 1^0,+) 
within the relevant timescale. The triplet states |ro,+) form a 
three-level system, whose dynamics is described by the mo- 
tion of a pseudo-spin- 1 Juj}. The conditional Hamiltonians in 



Eqs. ( [23] ) and ([24b]) are expressed in terms of the pseudo-spin 
as 



(m) 



(m) 



(25) 



The pseudo-spin model described in Eq. ( [25] ) is obtained in the 
zero magnetic field. The finite homogeneous magnetic field 
does not cause the mixing between singlet and triplet states. 



since it conserves the total spin. For the weak but non-zero 
magnetic field < y^^B < Dtj, the pseudo-spin model for 
the non-secular process is still valid but the average hyperfine 



field A|;^]j should be replaced by A'^j^ - y^B in Eq. (|25 ) 



With this pseudo-spin model, the electron spin decoherence 
(e.g., between | -h 1) and | - 1) electron spin states) induced by 
a nuclear spin pair under DD controls is expressed as 



L{ij}(t) 



1 1 



rTr 



(26) 

As an example, the nuclear spin pair induced decoherence un- 
der a 50-pulse PDD (PDD-50) control is shown in Fig. |5] 
Within the electron spin decoherence time < 6 ms un- 
der PDD-50 control, see the full calculation results for the 
decoherence time in Fig. [T3|, the pseudo-spin model repro- 
duces very well the exact result. The deviation from the ex- 
act result for larger time arises from the neglected dynamics 
driven by the hyperfine field diflference In particular, 

this pseudo-spin model shows that, in the weak field regime 
ynB < 1 Gauss, the secular spin flip-flop processes induced 
(24a) ^ component of a|^jj have negligible eflfect on the spin 

pair dynamics within the relevant timescale. 

Furthermore, the decoherence is sensitive to the dipolar in- 
teraction strength Dij. The dynamics of a separated nuclear 
spin pair [see Fig. [5jb)] is much slower than that of a dimer 
since the former has much weaker dipolar interaction. In Sec- 
tion [IV C 2[ we will show that the separated pairs have small 



contribution to the electron spin decoherence in the weak field 
regime. 

Thus, we conclude that, although both secular and non- 
secular processes can occur under weak magnetic field, the 
non-secular nuclear spin flipping is the dominant mechanism 
of electron spin decoherence, since it involves the average hy- 
perfine field of the two nuclear spins, while the secular process 
involves the much weaker hyperfine field diflTerence. 



2. Secular spin flip-flop 

When the Zeeman splitting is much larger than the dipolar 
interaction strength (y^B » Dij) and the hyperfine coupling 
(y^B » A^^^ and A^^^), the energy non-conserving process is 
suppressed by the Zeeman energy cost. The polarized states 
ITT) and are decoupled during the quantum evolution, and 
don't contribute to the electron spin decoherence (except for 
the static inhomogeneous broadening eflTect). The unpolar- 
ized states Itvl) and form a two-dimensional invariant sub- 
space [see Fig.|6];a)]. The nuclear spin pair in this subspace is 
mapped to a pseudo-spin- 1/2, i.e., Iti) i-^ IH) and ^ III)- 
Similar to the single nuclear spin case, the coherent evolution 
of the pseudo-spin is understood using the Bloch sphere pic- 
ture shown in Fig.jSJ^b). 

The nuclear spin flip-flop driven by the dipolar interaction 
gives rise to the transition between the pseudo-spin states Iff) 
and III). The transition rate is calculated through the nuclear 
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spin dipolar interaction as 



IV. RESULTS AND DISCUSSIONS 



Xij ^ (U l/^dipl U) = ^Dtj (l - 3 cos^ 6ij) , (27) 

where Oij is the angle between the pair orientation and the ex- 
ternal magnetic field. 

The diff'erence between the hyperfine fields projected along 
the direction of the external magnetic field induces an energy 
cost of the flip-flop 



Z^^ ^ (U \bn 



U) = z.(A|->-Ap. 



(28) 



Thus the flip-flop is mapped to the precession of a pseudo-spin 
(T about a pseudo-field h^"^^ = (x, 0, Z^'^A conditioned on the 
electrons spin state |m). 

The conditional Hamiltonian of the pseudo-spin reads 



im) 



(m) 



i(x,-cr, + ZyV,). 



(29) 



A typical evolution path of the pseudo-spin is shown in 
Fig.gb). 

Now that the pseudo-spin picture for the secular flip-flop 
of the nuclear spin pair is quite similar to the picture of sin- 
gle nuclear spin precession, the pair spin dynamics under DD 
controls can be understood in the similar way. The decoher- 
ence contributed by the nuclear spin pair in large magnetic 



field is obtained in the same form as in Eq. (20) with the 



eff'ective fields ^ replaced by the pseudo-fields h.'^s and 
the angle (fj replaced by the angle between the pseudo-fields 
(fij = arctan(Z,-^VX,7). 

In most cases, for the inter-nucleus distance of several 
angstroms, the flip-flop transition rate Xij (< 100 Hz) is much 
less than the energy cost Z^^^^ (~ kHz) due to the hyperfine 
field diff'erence. In this case, the pair contributes only a small 
fraction to the electron spin decoherence in the short time limit 
T ^ Xr.^ When the time gradually increases and the short 
time limit is not well satisfied, such kind of separated pairs 
will become important. A large number of such pairs give rise 
to a relatively smooth decay profile of the electron spin coher- 
ence. In Section IVC 3[ we will show that such kind of pair 
flip-flop process is the dominant decoherence mechanism in 
the strong magnetic field regime. 

In the opposite limit, i.e. Z^J^ <^ Xij, which corresponds 
to dimers located far away (e.g., several nanometers) from the 



NV center, two pseudo-fields h^^^ and h':J^^ are approximately 
parallel, i.e., the angle (ftj ^ 0. Such kind of remote dimers 
give small contribution to the electron spin coherence, and are 
not important under the strong magnetic field. 

Between these two limiting cases, a dimer not far 
(e.g., 1.5 nm) from the NV center can have dipolar interac- 
tion comparable to the hyperfine field diff'erence, i.e. Xtj ~ 
zf^^^ ~ kHz. Such dimers will induce coherent oscillations 
on the relatively smooth decoherence background. The ap- 
plication of the strong oscillation features in the atomic- scale 
magnetometry is discussed in Ref.Jli 



This section shows the calculation results of the NV cen- 
ter electron spin coherence. The NV center spin coherence 
has diff'erent decoherence behavior in diff'erent magnetic fields 
and at diff'erent timescales. The various decoherence behavior 
arises from the elementary processes discussed in the previ- 
ous Section. To identify the role of each elementary process 
in diff'erent magnetic fields, the following methods are used in 
analyzing the NV center spin coherence. 

First, in order to find out what kind of clusters are respon- 
sible to the electron spin decoherence, with the help of CCE, 
we resolve the contribution to the decoherence of nuclear spin 
clusters of diff'erent sizes. The CCE-1 will give the main con- 
tribution if the decoherence is caused mainly by the single 



spin dynamics described in Section III A While if the ^-spin 
correlation is the dominant decoherence mechanism, the CCE 
does not give the convergent result until the clusters of size n 
are included. Second, the role of diff'erent types of interaction, 
including the isotropic and anisotropic terms of the hyperfine 
coupling [see Eq. ([TO])], and the secular and non-secular terms 
of the nuclear spin dipolar interaction, are studied. To this end, 
we compare the exact results, which are obtained by solving 
the full Hamiltonian ([T]), with the results calculated by apply- 
ing the following approximations: (i) the isotropic approxi- 
mation, in which the anisotropic terms in hyperfine coupling 
such as SJ^. and S J^j are neglected, and (ii) the secular ap- 
proximation, in which the non-secular terms in nuclear spin 
dipolar interaction such as /^/| and /^/| are dropped. When 
both of the two approximations are applied, the only decoher- 
ence mechanism left is the secular spin flip-flop described in 



Section III B 2 Through the comparison, the importance of 
various types of interaction is identified. 



A. Free-induction decay 

In FID experiments, the coupling to the host ^^N nuclear 
spin induces fast oscillations (~ MHz) on the decoherence 
envelope induced by the ^^C bath spins. To single out the 
eff'ect due to the ^^C bath spins, only the envelope is presented 
here. 

Figure [TJa) shows the FID of an NV center spin in a typ- 
ical ^^C bath spin configuration under various magnetic field 
strength. The coherence decays in several yus. During such a 
short time, the dipolar interaction between nuclear spins (typ- 
ically < kHz) is too weak to take eff'ect. Thus, the FID is 
caused by the non-interacting nuclear spins, and the coherence 
is the product of the contributions of each single spins, i.e. 
^fid(0 = Y\ji-j,¥YD(t). This is confirmed in Fig. [TJb), where 
the CCE-1 results are convergent. Taking the pair correlation 
into account (i.e. CCE-2) does not change the coherence. 

The FID decoherence behavior depends on the magnetic 
field strength. The coherence time in strong field is longer 
than that in weak or zero field, and the decoherence profile 
deviates from the Gaussian decay in the medium field [e.g., 
B = 200 Gauss, see Fig.jTJa)]. 

It is the anisotropic nature of the hyperfine coupling in NV 
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FIG. 7. (a) FID of an NV center in a typical nuclear spin configura- 
tion under various magnetic field strength, (b) The contributions to 
FID of diff'erent cluster sizes (CCE orders), (c) The exact FID coher- 
ence and the coherence obtained by applying the isotropic approx- 
imation in zero and strong fields, (d) Histogram of FID coherence 
time T* distribution under zero and strong fields. 



centers that causes the distinguished decoherence features. In 
Fig. [TJc), we compare the exact FID coherence results with 
the results obtained with the isotropic approximation. In zero 
magnetic field, the anisotropic terms of the hyperfine coupling 
play an important role. If the anisotropic terms are ignored 
in calculations, the FID decoherence behavior in zero field 
will be identical to that in the strong field limit. While in the 
strong field case (e.g., B = 2000 Gauss), the contribution of 
the anisotropic terms is negligible. 

This decoherence behavior can be understood by analyzing 
the single nuclear spin dynamics. For each single spin, the 



contribution to the FID is shown in Eq. 20 ^a). In the weak 
field (B < I Gauss), the precession angle about the external 

a), the coherence is 



field ef ^ 0. According to Eq. 



20 



^fid(o ^ n 



cos ■ 



-r2?2/2 



(30) 



where A^^^^ is the magnitude of the hyperfine field of the 
7th nuclear spin with electron spin in | -h 1) state and F^ = 
2y (^^^^^) /4 is the FID decay rate in zero or weak fields. In 
the strong field limit (B » 300 Gauss), the eflTective fields h^^^ 

and h^.^^^ are almost parallel with each other, i.e. (fj ^ [see 
Fig.[3|. In this case, the coherence is 



M ^ n 



cos ■ 



of 



n 



cos 



(31) 



where A^.^^^ is the hyperfine field component along the mag- 



netic field direction, and F^ = 2j (^f/^) /4 is the FID deco- 
herence rate in strong filed limit. 



(+1) . 



always holds, the FID decoherence time 
in the strong field is longer than that in the weak field. Physi- 
cally, this is because, in the strong field limit, the ^^C nuclear 
spins are frozen along the quantization axis defined by the ap- 
plied magnetic field. Thermal fluctuations enable the nuclear 
spins pointing parallel or anti-parallel the magnetic field direc- 
tion, providing an Overhauser filed proportional to the hyper- 
fine fields components A^^^\ While in the zero field limit, the 
quantization axis of each nuclear spins can be chosen along 
the hyperfine field directions. The nuclear spins are parallel or 
anti-parallel along the hyperfine field, so that the Overhauser 
field strength is proportional to the magnitude A^.^^^ of the hy- 
perfine fields. The stronger Overhauser fields in weak mag- 
netic field induces the faster FID decoherence. In Fig. [7jd), 
we show the histogram of FID coherence time for about 200 
randomly generated bath spin configurations (the random po- 
sitions of 1.1% ^^C in the diamond lattice) under strong and 
zero magnetic fields. The FID decoherence time spreads over 
several microseconds for diflTerent spin configurations. The 
mean decoherence time over diflTerent configurations under 
zero and strong magnetic field is T* = V2/r = 2.1 yus and 
r*^ = V2/F2 = 3.6 yus, respectively. 

In a magnetic field between the strong and weak field limit 
(1 Gauss < B < 300 Gauss), as shown in Fig.jTJa), the FID 
decoherence profile deviates from the Gaussian decay. This 
decoherence behavior is caused by the quantum fluctuation 
induced by the ^^C spins close to the NV center. In a medium 
field, the Zeeman energy is comparable to the hyperfine cou- 
pling strength for the closely located ^^C spins. In this case, 
as shown in Fig.j3j the magnitudes of eflTective fields h^^^ and 

h^.^^^ are comparable and they are, in general, not parallel. 
Thus, the eigenstates of the noise operator Z?^ (i.e. the nuclear 
spin state has definite Overhauser field) are not eigenstates of 
the total Hamiltonian. In this way, the quantum fluctuation 
arises from the superposition of the nuclear spin states with 
diflTerent Overhauser fields, and causes the non-Gaussian de- 
cay in the medium magnetic field. This distinguished feature 
of FID decoherence is essentially attributed to the anisotropic 
nature of the hyperfine coupling of NV centers. In QDs and 
shallow donors, where the isotropic hyperfine coupling dom- 
inates, the external field and the Overhauser field are always 
parallel, and such non-Gaussian FID would be absent. 



B. Hahn echo 

The thermal noise is eliminated by applying the Hahn echo. 
The decay of the spin echo is caused by the dynamical fluc- 
tuations of the nuclear spin bath. The decoherence time of a 
certain NV center spin depends on the bath spin configurations 
around it.^^ In Fig.]8[a), we show a typical Hahn echo coher- 
ence as a function of the total evolution time and the magnetic 
field strength. 

Three magnetic field regimes, in which the NV center spin 
coherence decays in diflTerent manners, are roughly distin- 
guished by the horizontal lines in Fig.[8ja). In the weak field 
regime (B < I Gauss), the coherence decays monotonically 
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FIG. 8. (a) Hahn echo of NV center spin coherence between |0> and | + 1) as a function of total evolution time t and magnetic field B. 
The strong, medium and weak magnetic field regimes, in which the NV center spin has diff'erent decoherence behavior, are separated by the 
horizontal lines. The solid (dashed) line indicates the time at which the NV center spin coherence decays to the Xje value for the first peak 
(envelope), i.e. the time Tw (^2) defined in (b). (b)-(f) Cross-sections of (a) for various magnetic field strength, (b) The collapse and revival 
eff'ect under 30 Gauss magnetic field with three diff'erent timescales indicated. The red dashed line is the envelope of the revival peaks, which 
decays to 1/^ at time T2. Inset: close-up of the first several peaks. The half peak width Tw and the revival period Tr are indicated, (c) In a 
weak magnetic field of 0.3 Gauss (approximately the earth magnetic field strength), the NV center spin coherence decays in a timescale Tw, 
and does not revive, (d)-(f) The suppression of the collapse and revival eff'ect by increasing the magnetic field strength. For 5 = 100 Gauss 
in (d), the revival peaks get closer to each other. The coherence does not collapse to zero and the revival peaks are not well resolved under 
B - 500 Gauss in (e). For B - 3000 Gauss in (f), the collapse and revival eff'ect is completely suppressed. 



within ~ 200 yus [see Fig. [SJc)]. In the medium field regime 
(1 Gauss < B ^ 300 Gauss), the main feature is the peri- 
odic coherence collapse and revival under an overall decay 
envelope [see Fig. [8jb) and (d)]. In the strong field regime 
(B > 300 Gauss), the collapse and revival effect is greatly 
suppressed [Fig. [8je)] and finally vanished [Fig. [Sjf)]. The 
coherence monotonically decays as in the weak field regime, 
but with a longer timescale of ~ 700 yus. 

To quantitatively characterize the decoherence behavior, 
three timescales are defined [see Fig. [8jb)]: (i) the electron 
spin coherence collapses and revives within a period Tr, form- 
ing a periodic peak structure, (ii) for each coherence peak, the 
half width is denoted as Tw (the half width at 1/^ height, ap- 
proximately the same for all the peaks), and (iii) the overall 
envelope decays to 1 in a slower timescale T2 ~ ms. The 
peak decay time Tw of the first (half) peak and the envelope 
decay time T2 are highlighted in Fig. [Sja) with the solid and 
dashed contours, respectively. 

The field dependence of Tr, Tw and T2 are summarized in 
Fig. [9] The revival period Tr is inversely proportional to the 
magnetic field strength as expected. The envelope decay time 
T2 increases with increasing the magnetic field in the medium 
field regime, and becomes constant in the strong field regime. 
The peak decay time Tw decreases with decreasing the mag- 
netic field as Tw oc B~^ with k ^ 0.6 in the medium field 



regime, and approaches the revival period Tr when the field is 
close to the boundary of the strong field regime. In the weak 
field regime, deviates from the B~^ law, and is saturated 
in the zero field limit. The magnetic field regimes are divided 
according to not only the apparent field dependence of Tr, 
Tw and T2 shown in Fig. [9j but also the different underlying 
decoherence mechanisms, which will be discussed below. 



7. Medium magnetic field regime 

Let us start from the medium field regime with the mag- 
netic field 1 Gauss < B < 300 Gauss. The contribution to de- 



coherence of different cluster sizes is shown in Fig. 10 c). As 
discussed in Section |lll| the collapse of the electron spin co- 
herence is induced by the single ^^C nuclear spin precession. 



This is confirmed in Fig. 10c), where the CCE-1 gives rise to 
the periodically collapse and revival. The coherence would re- 
cover perfectly if only the single-spin clusters are considered. 
The envelope decays when the two- spin clusters are taken into 
account (CCE-2). This means the envelope decay is caused by 
nuclear spin pair dynamics. The higher order correlations of 
more than two nuclear spins are not important in Hahn echo 
since the CCE-2 gives the convergent result. 

Furthermore, we identify the contributions to decoherence 
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FIG. 9. The magnetic field dependence of Tw, and The black 
line with square symbols is the calculated peak width, and the red 
solid line is obtained by fitting the data in the region 1 Gauss < B < 
100 Gauss. The peak width Tw depends on B as B~^ with k - 0.63. 
The value of k slightly depends on the specific nuclear spin positions 
around the NV center. For B < \ Gauss, the peak width devi- 
ates from the B~^ law, which defines the weak magnetic field regime. 
The revival time Tr is inversely proportional to the magnetic field 
strength B (circle symbols and blue dashed line). The rough bound- 
ary between the medium and strong magnetic field regimes is defined 
by the approaching of Tr to about 5 = 300 Gauss. The envelope 
decay time (black line with diamond symbols) in the strong mag- 
netic field (> 300 Gauss) is constant, and decreases as the magnetic 
field is decreased to ~ 10 Gauss. For 5 < 10 Gauss, there are too few 
revival peaks to determine the envelope decay time, and Ty^ should 
be taken as the relevant decoherence time. 



of different types of interaction. As shown in Fig. 10 ^d), the 
secular approximation gives almost the same result as the 
exact one. This implies the non-secular spin flipping has 
negligible effect on spin decoherence in the medium field 
regime. It is understandable since the dipolar interaction 
strength (<^ kHz) is much smaller than the Zeeman energy 
in this field regime (> kHz). Since the anisotropic coupling 
of the hyperfine interaction is responsible to the collapse and 
revival eflfect as discussed in Section IHTaI the coherence col- 
lapse and revival do not occur when the isotropic approxima- 
tion is applied. In addition, the anisotropic terms results in a 
shorter envelope decay time than the coherence time with 
only the secular flip-flop process taking into account [i.e. with 
both secular and isotropic approximation applied, see the blue 
dashed line Fig.[TO];d)]. 

The effect of the anisotropic terms on the coherence enve- 
lope is explained as follows. Notice that, besides the secular 
flip-flop process | tvl) ^ I vlT). the energy non-conserving 
two-spin flipping processes, such as | TT) ^ I \X)-> are induced 
by the transverse hyperfine field Aj^^ or ^j,y In the medium 
field regime, these energy non-conserving processes can oc- 
cur since the hyperfine interaction strength for the nuclear pair 
close to the NV center (> 10 kHz) could be stronger than or 
comparable to the nuclear spin Zeeman energy. With such 
kind of processes activated, the envelope decay in the medium 
field regime is faster than the decoherence solely caused by 
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FIG. 10. Decoherence mechanisms of Hahn echo in different mag- 
netic field regimes, (a) The contributions to decoherence of different 
cluster sizes (CCE orders) in 5 = 0.3 Gauss, (b) The contributions 
to decoherence of different types of interaction m B - 0.3 Gauss. 
The red line with square symbols is the coherence calculated with 
the isotropic approximation of the hyperfine interaction. The green 
line with diamond symbols is the coherence calculated with the sec- 
ular approximation. The blue dashed line is the coherence obtained 
with both isotropic and secular approximations. The black solid line 
indicates the exact results without approximations, (c) and (d) are 
the same as (a) and (b), respectively, but for 5 = 15 Gauss, (e) and 
(f) are the same as (a) and (b), respectively, but for B - 3000 Gauss. 



the secular nuclear spin flip-flop process [the blue line in 
Fig.[TOjd)]. With increasing the magnetic field strength, this 
kind of processes are gradually suppressed due to the increas- 
ing Zeeman energy. As a result, the envelope decay time 
[the dashed contour in Fig. [8ja) and the diamond symbols in 
Fig. |9j increases and approaches the value in the strong field 
regime. 

As for the single peak decay time Tw, we plot the magnetic 
field dependence of Tw in Fig.|9]in comparison with the re- 
vival period Tr. According to Eq. ( [2Q| ), in the short time limit, 
the electron spin coherence is expressed as 



(32) 



which implies the decay time of a single peak scales with the 
magnetic field strength as ~ 5"^/^. The numerical result gives 
the peak width Tw ^ with the index k ^ 0.6. The devi- 
ation of index k from 1 / 2 is due to the violation of the short- 
time condition for the ^^C spins close to the NV center. 
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2. Weak magnetic field regime 

The revival period Tr increases as the magnetic field is de- 
creased. As shown in Fig. [8] for 5 < 2 Gauss, the subsequent 
revival peaks are hardly observed due to the decay of the en- 
velope, and only the first (half) peak is relevant in this case. 
Meanwhile, for 5 < 1 Gauss, the dependence of Tw on B be- 
gins to deviate from the law as governed in the medium 
magnetic field regime, and becomes saturated approaching the 
zero field (see Fig. [9]). Thus, B = I Gauss defines the rough 
boundary of the weak field regime. In this regime, the deco- 
herence time is defined by Tw 

To analyze the decoherence mechanism in the weak field 



regime, in Fig. 10 ^a), we check the CCE convergence to find 
out the relevant clusters responsible to the decoherence. Being 
diff'erent from the medium field regime, where the single spin 
dynamics dominates the peak decay, the single spin preces- 
sion process becomes less important in the weak field regime 



(see Section III A). As shown in Fig. 10 a), under a mag- 
netic field B = 03 Gauss, the single spin precession (CCE-1) 
contributes a part of the decoherence, and the coherence con- 
verges when the pair dynamics (CCE-2) is taken into account. 

In Fig.[TOjb), diff'erent types of interaction in the weak field 
regime are analyzed. In this regime, the conservation of Zee- 
man energy is not a necessary requirement. Thus, the Zee- 
man energy non-conserving spin flipping processes induced 
by anisotropic hyperfine coupling and the non- secular terms 
in the dipolar interaction are activated. As discussed in Sec- 
tion |IIl| the noise amplitudes produced by these processes are 



typically stronger than that produced by the secular flip-flop 
process of nuclear spin pair. In this way, the energy non- 
conserving spin flipping processes become the dominant de- 
coherence mechanism [see the red and blue lines with symbols 
in Fig.[TOjb)], while the secular flip-flop process has negligi- 
ble contribution to the decoherence [see the blue dashed line 
inFig.lTgb)]. 



3. Strong magnetic field regime 

Since Tw decreases slower than Tr, as shown in Fig. |9] in 
strong magnetic fields {B > 300 Gauss), Tr and Tw are of 
the same order. The collapse and revival eff'ect due to the sin- 
gle nuclear spin precession is greatly suppressed in the strong 
field regime, and the peak structure will finally disappear in 
the strong field limit [see Fig.[8jd)-(f)]. 

In Fig. 10 e), the CCE convergence is shown. In strong 
magnetic field, the nuclear spins are frozen along the mag- 
netic field direction, and the single spin precession (CCE-1) is 
almost completely suppressed. The decoherence is caused by 
the nuclear spin pair flip-flops (CCE-2). Furthermore, all the 
energy non-conserving two-spin processes are suppressed by 
the strong field, leaving the secular flip-flop process being the 
only decoherence mechanism in the strong field regime. This 
is confirmed in Fig. 10 ^f), where neglecting the anisotropic 
coupling and the non- secular terms does not aff'ect the de- 
coherence behavior. Since the absolute magnitude of Zee- 
man energy is not involved in the secular flip-flop process. 



the coherence time T2 does not depends on the magnetic field 
strength anymore in the strong field regime [see the dashed 
contour in Fig. [5] and the diamond symbols in Fig.[9|. 



C. Multi-pulse dynamical decoupling controls 

In this Section, we study the eff'ects of diff'erent decoher- 
ence mechanisms in diff'erent magnetic field regimes, for NV 
center spins under multi-pulse DD control. 



1. Medium magnetic field regime 

Let us first consider the medium magnetic field regime, in 
which the spin coherence under DD controls has been studied 
experimentally.^^ As a direct generalization of the Hahn echo, 
the electron spin coherence under PDD control also presents 
the collapse and revival effect [see Fig.[TTJa)]. We will first 
focus on the dynamics within a revival peak, and investigate 
the DD control eff'ect on the peak decay time. 

The electron spin coherence under magnetic field B = 
60 Gauss and multi-pulse PDD controls is shown in Fig.[TTJb). 
The peak decay time Tw is extended from ~ 10 yus for the 
Hahn echo to ~ 600 yus for the 100-pulse PDD control. The 
small oscillations on the decoherence curves arise from some 
^^C nuclear spins located relatively close to the NV center 
which have strong hyperfine interaction (A^.^^^ ~ 100 kHz). 
The coherence time increases linearly with increasing the DD 
control pulse number [see the inset of Fig. [TTJb)], which 
is consistent with the experimental observations in the same 
magnetic field regime.^ In contrast, the experiment of NV 
center coherence in electron spin baths shows that the co- 
herence time depends on the DD control pulse number as 
^2/3 23 jj^jg |g because the diff'erent features of the noise spec- 
tra in nuclear and electron spin baths. In electron spin baths, 
the noise spectrum is of Lorentzian shape with a long high- 
frequency tail, while the noise spectrum due to the nuclear 
spin bath has a hard high-frequency cut-off'. In particular, in 
the medium magnetic field regime, the noise spectrum is a 
peak located around the Lammor frequency. The high fre- 
quency noises in electron spin baths make the DD less efficient 
in protecting the NV center spin coherence (i.e. the coherence 
time increases slower with increasing pulse number). The 
quantitative relation between the coherence and noise spec- 
trum is an interesting problem but beyond the scope of this 
paper. It will be discussed elsewhere. 

As shown in Fig. [TTJc), the peak decay time under multi- 
pulse PDD controls has the similar magnetic ffeld dependence 
as in the Hahn echo case. The similar magnetic ffeld depen- 
dence of the peak width suggests that the same decoherence 
mechanism, i.e. the single spin precession process, is respon- 
sible for the single peak decay under multi-pulse PDD con- 
trols. To conffrm this point, the CCE convergence for the 
100-pulse PDD control is shown in Fig.[TTJd). It is observed 
that CCE-1 almost gives the convergent coherence. Including 
the correlations between more spins induces negligible change 
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FIG. 11. (a) Electron spin coherence between |0> and | + 1> un- 
der PDD-10 control in the medium magnetic field regime (B = 
60 Gauss). Similar collapse and revival eff'ect is observed, but with 
a longer timescale than the Hahn echo case, (b) The decoherence 
within the first (half) peak under PDD controls of various pulse num- 
bers indicated by the integers associated with each curves. For the 
sake of clarity, the subsequent revivals are not shown. The inset 
shows that the decay time, at which the coherence decays to l/e 
value, increases linearly with increasing pulse number, (c) The mag- 
netic field dependence of the decay time (at l/e height) for various 
PDD pulse numbers, (d) The CCE convergence of the PDD- 100 re- 
sult under B = 60 Gauss. The single spin dynamics (CCE-1) is the 
dominating decoherence mechanism in this magnetic field regime. 



to the coherence result. This convergence indicates that, al- 
though the coherence time under multi-pulse DD is prolonged 
by more than one order of magnitude as compared with the 
Hahn echo, the dominating mechanism for the initial coher- 
ence decay is still the single spin precession. 
In Fig. 



T2[a) and (b), we check the CCE convergence in a 
longer timescale to study the envelope decay. Similar to the 
Hahn echo case, the single-spin precession (CCE-1) results 
in the perfect revival, while the nuclear spin pair correlation 
(CCE-2) contributes the main part of the envelope decay. As 
the coherence is prolonged under DD control, the higher or- 
der correlations (clusters with more than two spins) begin to 
take effect on the envelope decay and give a small correction 



to the CCE-2 result. As shown in Fig. 12 b), the calculation 
converges with the correlations of six- spin clusters included 
(CCE-6). 
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FIG. 12. (a) and (b) The contributions to overall decoherence of 
diff'erent cluster sizes (CCE orders) in 5 = 20 Gauss and PDD-10 
control. CCE-1 give rise to the collapse and revival eff'ect. The enve- 
lope decay is caused mainly by CCE-2. The calculation is converged 
at CCE-6. (c) The contributions to decoherence of diff'erent types of 
interaction in 5 = 20 Gauss and PDD-10 control. 



The contribution of diflTerent type of interaction under DD 



control is shown in Fig. 12 c). The conclusion is similar to 
the Hahn echo case. In the medium field regime, the non- 
secular spin flipping contributes little to the decoherence. The 
anisotropic terms of hyperfine coupling cause the collapse and 
revival phenomenon and provide the additional envelope de- 
cay besides the one caused by the secular spin flip-flop pro- 
cess. 



2. Weak magnetic field regime 

In the weak magnetic field regime {B < \ Gauss), the two 
electron spin states | ± 1) are hardly resolved through the tran- 
sition frequencies. In this case, instead of the coherence be- 
tween |0) and 1 + 1) states, we investigate the coherence be- 
tween 1 + 1) and I - 1), which can be generated and manipulated 
in experiment by linearly polarized microwave pulses. ^^^^ 

The electron spin coherence under 5 = 0.3 Gauss and 
multi-pulse PDD controls is presented in Fig. 



13 Similar to 



the medium field case, the coherence time Tw increases lin- 



early with increasing pulse number (see the inset of Fig. 13 ). 
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FIG. 13. Electron spin coherence between | + 1) and | - 1> under weak 
magnetic field (0.3 Gauss) and PDD controls. The integer numbers 
associated with each curves are the pulse numbers of the PDD con- 
trols. Inset: the coherence time (at 1/^ height), increases linearly 
with increasing pulse number. 



To analyze the decoherence mechanisms, we examine the 
CCE convergence. Taking the coherence under PDD-50 con- 



trol for example [Fig. 14 a)], the CCE-1 has a small contri- 
bution to the decoherence, and the CCE-2 almost produces 
the converged results. This means the dominant decoherence 
mechanism under weak magnetic field is the nuclear pair cor- 
relations. Higher order correlations of larger nuclear spin 
clusters, e.g., the CCE-5 and CCE-6, provide only small cor- 
rections in details. 

We further categorize all the nuclear spin pairs into two 
classes according to the inter-nuclei distance. The first kind is 
the nuclear spin dimers. For natural abundance ^^C, there are 
on average about 12 dimers in the bath within 4 nm around the 
NV center. The second class of pairs consists of a large num- 
ber of separated spin pairs, which have the inter- spin separa- 



tions larger than the C-C bond length. In Fig. 14b), the diflfer- 
ent eflTects of these two classes of pairs are shown. The deco- 
herence in the weak field regime is dominated by the nuclear 
dimers [see Fig.|5ja) for decoherence by a typical dimer]. The 
separated pairs, which have much slower dynamics than the 
dimers because of the weaker dipolar interaction, only give a 
small contribution to the total decoherence [see Fig. |5jb) for 
decoherence induced by a typical separated pair] . 



In Fig. 14 c), the roles of different types of interaction 
under DD control are compared. Similar to the Hahn echo 
case, the anisotropic terms of hyperfine coupling and the non- 
secular part of the nuclear spin dipolar interaction are the rele- 
vant decoherence mechanisms under DD controls in the weak 
field regime. The secular flip-flop process contributes only 
a small amount of decoherence [see the blue dashed line in 
Fig. 14 c)]. This is consistent with the discussion in Sec- 
tion llllBll 
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FIG. 14. (a) The CCE convergence for PDD-50 under weak magnetic 
field B = 0.3 Gauss. The 2-spin clusters give the dominating contri- 
bution (CCE-2), and the result converges by including up to 6- spin 
clusters (CCE-6). (b) Contributions to NV center spin decoherence 
from ~ 12 nuclear spin dimers (red solid line) and the large number 
of separated pairs (blue dashed line) in the bath. The decoherence in 
the weak field regime is mainly caused by the dimers. (c) Decoher- 
ence caused by diff'erent types of interaction in B = 0.3 Gauss and 
PDD-50 control. 



3. Strong magnetic field regime 

Finally, we study the decoherence in the strong magnetic 
field regime. The coherence under multi-pulse PDD controls 



is shown in Fig 15 ^a). The overall coherence time T2 is lin- 
early increased with increasing the pulse number [see the in- 



set of Fig 15 a)]. Meanwhile, in the same strong field (e.g., 
B = 3000 Gauss), sharp dips appear on the smooth back- 
ground decay, and the dips get more profound as increasing 



the control pulse number [see Fig.[T5[b)]. 

These coherence dips under higher order DD control come 
from the noises produced by the single spin precession. Al- 
though most of the nuclear spins are frozen by the large Zee- 
man energy in the strong magnetic field, a few ^^C spins lo- 
cated close to the NV center (with the hyperfine coupling 
~ 100 kHz) could produce weak noises to the NV center spin. 
In the Hahn echo case, such weak noises can hardly cause any 
notable effect on the coherence. While in the multi-pulse PDD 
cases, the weak noises can be greatly amplified by the pulse 
sequence, and cause profound coherence dips. Under PDD 
controls, the coherence dips are quasi-periodic. The period is 
proportional to the pulse number N, and the dip depth is pro- 
portional to [see Fig. 



15 



^b) and Ref. 11 1. This weak noise 
amplification by multi-pulse DD control has the potential ap- 
plication in the single nuclear spin detection. The detailed 
discussion on this prospect will be presented elsewhere. 

The origin of the coherence dips is confirmed by the CCE 
convergence check. As shown in Fig[T6ta), under PDD-4 con- 



15 



CD 
O 

c 

CD 
O 

O 



CD 
O 

c 

CD 



O 
O 




0.00 



0.01 



0.02 



0.03 



Time (ms) 



FIG. 15. (a) Electron spin coherence between |0> and | + 1> under 
strong magnetic field B = 3000 Gauss and PDD controls. The in- 
teger numbers associated with each curves are the pulse numbers of 
the DD controls. Inset: the coherence time T2 linearly increases as 
increasing the pulse number, (b) The close-up of the initial stage of 
the decoherence under PDD-5 and PDD- 10 controls. Sharp coher- 
ence dips are caused by the noise from the nuclear spins close to NV 
center. 
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trol, the single spin dynamics (CCE-1) does not cause the co- 
herence decay but provides the coherence dips. Similar to the 
weak field regime, the decoherence is mainly caused by the 
pair correlations (CCE-2), while the higher order spin corre- 
lations induce corrections to the details. 

Furthermore, we compare the contributions of the dimers 
and the separated pairs. As we have shown in our previous 
study,^ a dimer located close enough to the NV center (with 
distance < 1.5 nm) can induce strong oscillations of the elec- 
tron spin coherence. In this paper, we are mainly interested in 
the overall decoherence, which is contributed by many weakly 
coupled pairs in the bath. As a result, in Fig. [T6j we choose 
a bath spin configure without a strongly coupled dimer within 
1.5 nm to NV center (about 50% probability to find such kind 
of bath spin configurations in diamond with natural isotope 
abundance) to get relative smooth background decoherence 
profile. As shown in Fig [T6jb), in contrast to the weak field 
regime, the decoherence under strong magnetic field is domi- 
nated by the separated pairs. 

The diflTerent roles of dimers in decoherence under weak 
and strong magnetic fields can be understood as follows. Be- 
cause of the energy conservation, only the secular spin flip- 
flop process is eflTective in strong magnetic field. As discussed 



in Section |IIIB[ the secular spin flip-flop process is associ- 
ated with the hyperfine field diflTerence between the two nu- 
clear spins. Because of the small hyperfine field diflTerences 
for dimers far from the NV center, the two pseudo-fields h^^^ 
are approximately parallel with each other. Thus, according 
to the pseudo-spin picture in Fig.[6j the remote dimers are not 



FIG. 16. (a) The CCE convergence for PDD-4 under strong magnetic 
field B = 3000 Gauss. The 2- spin clusters give the main contribution 
(CCE-2), and the result converges by including up to 6-spin clusters 
(CCE-6). Inset: close-up for the CCE-1 and CCE-2 contributions 
to the coherence dips in short time, (b) Contributions to the deco- 
herence from ~ 12 nuclear spin dimers (red solid line) and the large 
number of separated pairs (blue dashed line) in the bath. The electron 
spin decoherence under strong magnetic field is mainly caused by the 
separated pairs, (c) Decoherence caused by diff'erent types of inter- 
action under B = 3000 Gauss and PDD-4 control. Inset: close-up for 
the contributions to the coherence dips of diff'erent interactions. 



important for the decoherence. This discussion is supported 
by the analysis of the contributions of different type of inter- 



action shown in Fig. 16 c). Except the coherence dips induced 



by the anisotropic terms [see the inset of Fig.[T6[c)], neglect- 
ing the energy non-conserving spin flipping processes does 
not change the NV center spin decoherence behavior. Thus, 
we conclude that the secular flip-flop of the separated nuclear 
spin pairs is the main decoherence mechanism in the strong 
magnetic field regime. 



V. CONCLUSION 

In this paper, we systematically investigate the electron spin 
decoherence of NV centers in ^^C nuclear spin baths under 
DD controls. We present the different decoherence behaviors 
in weak, medium and strong magnetic field regimes, and ex- 
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amine the underlying microscopic mechanisms. In medium 
magnetic field regime, the single nuclear spin precession pro- 
cess induces the coherence collapse and revival eff'ect. In 
the weak and strong magnetic field regimes, the decoherence 
is dominated by two kinds of nuclear spin pair processes, 
namely, the non-secular spin flipping of dimers in the weak 
field regime and the secular spin flip-flops of the separated 
nuclear spin pairs in the strong field regime. 

The NV center systems investigated here are diff'erent from 
QDs and shallow donors by the dipolar nature of the hyperfine 
interaction and the diluteness of the nuclear spin bath. The NV 
center electron spin in the nuclear spin bath is a representative 
system for deep-level defects. The decoherence behavior of 
NV centers discussed here can be generalized to other sim- 
ilar systems."^^ Better understanding of the underlying deco- 
herence mechanisms will benefit various applications of such 
kind of systems, such as quantum information processing!^ 
and ultra- sensitive metrology."^^ 




FIG. 17. Geometric picture of the effective bath spin precession 



about the axes along n^~^ with precession angles 6)[\ under PDD- 
^ > 1) control. The angle between the precession axes is lak. The 
vectors sin(0^~V2) are decomposed into the common part and the 
difference part r^^. 



- o.^i^) withn^*^ 



where the effective precession angle 6]^^ = 6ih^ vvi..xx 
being the directions of the effective precession axis and the 
magnitude 6i is determined b}E51 
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The vector of effective precession angles are decomposed into 
the common part Ri and the difference part ri as 



(±) 



sin = Ri ± ri. 



(A.5) 



Appendix: Single-spin dynamics under FDD controls 

Here we present the derivation of the electron spin deco- 
herence induced by single spin dynamics under FDD con- 
trols. To be specific, we will use the terminology for single- 
spin dynamics, but the discussion automatically applies to the 
pseudo-spin dynamics. Given the electron spin in state |m), 
the free evolution operator of a nuclear spin is denoted as 



-10': 



•^/^ = cos 



nim) 
A 

2 
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-/tr-sin-^, (A.l) 



□ (m) 



where sin(v) = (v/v) sin(v), and the angle O""^^ = y^M^^T is 
the precession an gle f or the nuclear spin about the effective 
field h^'^^ [see Eq. (18)]. Here we use the integer subscript k to 
indicate the DD control order = means free evolution, and 
1 for Hahn echo), while the nuclear spin index j is omitted for 



the sake of clarity. According to Eq. ( 15 ), the single nuclear 
spin contribution to Hahn echo is 



^Hahn (2t) = ^Tr 



{U[^\r)f u[-\r) 



(A.2) 



Here we have defined the effective evolution operators 
under the Hahn echo control as 
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(A.6b) 



The geometric picture of the effective precession is shown in 



Fig. ([17]). Notice that the common part Ri and the difference 
part ri are perpendicular to each other. The angle between the 
two effective precession axes is determined by 



tanofi 



IRil 



(A.7) 



substituting Eqs. (|A.3|)-( A.6) into Eq. (A.2), one obtains 



LHahn(2T)= l-2|riP 



2 ^0 



= 1-2 sin^ (fi sin^ — sm 

^2 2 



(A.8) 



where (p is the angle between the effective fields h^^^ and h^^^\ 
The effective evolution operators for higher order FDD 
controls are defined by combining the evolution operators of 
lower orders in the similar way. For the even order FDD con- 
trol FDD-2^, the effective evolution operators U^\t) are 



U'^\T) = [uf\T)Uf(T)f^\e 



-16^ 



cr/2 



17 



where O^^j^ = 02k^2k ~ ^^2 ~ ^^2^2 • precession angle 
62 can be obtained through 
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2 2 2 
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(A. 10a) 
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The directions of precession axes are determined by the vec- 
tor sin(^2^V2), which are decomposed into the common and 
difference parts in the similar way of Eq. \A.5) as 



e 



(±) 

2k 



sin =^2k± Y2k' 



For k= 1, the vectors R2 and Y2 are calculated as 



n 

R2 = 2cos— Ri, 
2 

Y2 = 2Ri X ri. 



(A. 11) 



(A. 12a) 
(A. 12b) 



For k > 1, the vectors sin(^2lV2) have the same direction 

as sin(^2~V2), but with a magnitude sin(^^2/2). The angle 
between effective precession axes is determined by 



tanQf2 = 



|ril 



IR2 



cos 61/2 



(A.13) 



which is the result shown in Eq. (22a). 



With Eqs. ( |A.9| )-( |A.13) , the coherence under PDD-2^ con- 
trol is expressed in terms of the angles 62 and 0^2 as 



L2k('^kT) = -Tr 



= 1 - 2\Y2k\ =1-2 sin^ a2 sin^ 



ke2 



(A. 14a) 
(A. 14b) 



The odd order PDD control PDD-(2^ -h 1) is a combina- 
tion of PDD-2^ and a Hahn echo. The evolution operator is 
^S+i = uf\T)Ufj^{T\ and the coherence at t = {Ak + 2)t is 

L2k.i{t) = ^Tr [(t/^;^i(T))' U^Ut)] (A.15a) 

= \Tx[Q\{T)Q2k{T)\, (A.15b) 

where 

21 (r) = [U^;\T)j U^^\t) = LHahn(2T) - iCT • ^1, 

(A. 16a) 

Q2k{T) = U^-,\t) {U^2k(r)f = L2k(4kT) - ia • 6^, 

(A. 16b) 

with the traceless terms determined by the vectors 



^1=2 (vi X Ri cos y ri j , (A. 17a) 

6^ = 2 (r2;^ X R2;^ - cos ^r2k) . (A.17b) 



Then, the coherence L2k+i(t) is the sum of LHahn(2T)L2y^(4^T) 
and the correction term 



^corr(T) - ' ^2k- 



(A.18) 
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